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O ■ Abstract 

We propose a general framework that leads to one-dimensional XX and Hubbard 
^ ' models in full generality, based on the decomposition of an arbitrary vector space (pos- 

sibly infinite dimensional) into a direct sum of two subspaces, the two corresponding 
orthogonal projectors allowing one to define a i?-matrix of a universal XX model, and 
then of a Hubbard model using a Shastry type construction. The QISM approach 
ensures integr ability of the models, the properties of the obtained iZ-matrices leading 
to local Hubbard-like Hamiltonians. 

In all cases, the energies, the symmetry algebras and the scattering matrices are 
explicitly determined. The computation of the Bethe Ansatz equations for some sub- 
sectors of the universal Hubbard theories are determined, while they are fully computed 
in the XX case. A perturbative calculation in the large coupling regime is also done 
for the universal Hubbard models. 
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1 Introduction 



The celebrated Hubbard model, introduced in the sixties [1,2] in order to study strongly 
correlated electrons, has been widely studied since then, essentially in connection with con- 
densed matter physics. Due to the extent of the literature on the subject, the reader is 
invited to refer to the books [3,4] and references therein. The eigenf unctions and energies 
of the ID-model are known by means of the Bethe ansatz thanks to the works of Lieb and 
Wu [5], the complete set of eigenstates being obtained in [6], exploiting the SO (4) symmetry 
present in the one-dimensional case. 

The essence of the Hubbard model is rather fascinating: although the one-dimensional 
model was solved in the late sixties, the understanding of the model in the light of the 
quantum inverse scattering method became clear only twenty years after with the works 
of Shastry [7,8] and Olmedilla et al. [9]. The main idea is to couple the R-matrices of two 
independent XX models, through a term depending on the coupling constant of the Hubbard 
potential. The complete proof of the Yang-Baxter relation for the Hubbard R-matrix was 
given by Shiroishi and Wadati [10]. 

Since then, generalizations of the Hubbard model in the framework of the i?-matrix 
formalism have been proposed. A first step was done by Maassarani [11, 12], extending the 
i?-matrix construction to the gl{N) case. 

The appearance of the Hubbard model in the context of N = 4 super Yang-Mills the- 
ory led to new motivations to investigate further the supercase. The Hubbard model at 
half-filling, when treated perturbatively in the coupling [13], reproduces the long-ranged in- 
tegrable spin chain of Ref. [14] as an effective theory. It was conjectured in [14] that the 
Hamiltonian of this chain be an all-order description of the dilatation operator of iV = 4 
super Yang-Mills in the su(2) subsector. There may be the possibility that some integrable 
extension of the Hubbard model (e.g. involving superalgebras) could be put in relation to 
other subsectors of the N = 4 super Yang-Mills theory. 

The Hubbard model has also arisen in [15], where an S'-matrix for a long-range interacting 
integrable quantum spin chain with centrally extended sii(2|2) symmetry was constructed. 
This S'-matrix was shown to be proportional to Shastry's i?-matrix up to a dressing phase. 
This phase indeed leads to a breakdown of the conjecture of [14] beyond three loops and to 
transcendantal contributions to the dilatation operator eigenvalues. However, the proposal of 
a string Bethe ansatz and the appearance of the Hubbard i?-matrix in the study of integrable 
structures in view of the AdS / CFT correspondence ask for learning more about generalized 
Hubbard models. Hopefully, the statistical mechanics community may also find interest in 
exploring these structures. 

A superalgebraic generalization of the Hubbard model in the spirit of Shastry's con- 
struction has been proposed by the authors in [16] and [17], where a general approach to 
constructing a number of super Hubbard models was developed. Each of the obtained mod- 
els can be treated perturbatively and thus gives rise to an integrable long-ranged spin chain 
as an effective theory. The symmetry of the super Hubbard model based on gl(m\n) was 
shown to be gl(m — l|n — 1) © <7^(1|1) © gl( m — l|n — 1) © gl(l\l). In this paper, we propose 
a general framework that leads to XX and Hubbard models in full generality. It may also 
constitute an interesting starting point for dealing with integrable bosonic Hubbard mod- 
els. More precisely, it is based on the decomposition of an arbitrary vector space (possibly 
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infinite dimensional) into a direct sum of two subspaces, the two corresponding orthogonal 
projectors allowing one to define a i?-matrix of a universal XX model, and then of a Hubbard 
model using a Shastry type construction. The QISM approach ensures the integrability of 
the models, the properties of the obtained i?-matrices leading to local Hubbard-like Hamil- 
tonians. In the finite dimensional case, they can be interpreted in terms of 'electrons' after 
a Jordan- Wigner transformation [18] (see some examples in [16]). 

The plan of the paper is as follows. In section [2j we extend the construction of XX 
models for algebras [19] and superalgebras [16] to the case of an arbitrary vector space, 
possibly infinite dimensional. We focus to the general case g/(m|n) in section [31 in which 
the Hamiltonians are explicitly constructed and the Bethe ansatz equations computed. In 
section HI we tackle with the case of universal Hubbard models, performing in the same 
way the calculation of the i?-matrices, the transfer matrices and corresponding Hamiltoni- 
ans. For both models (XX and Hubbard), the energies and the symmetry algebra, which 
is related to the choice of the projectors, are determined, and the corresponding charges 
computed. Section is devoted to the Bethe ansatz equations for universal gl(vn\n) Hub- 
bard models. The computation of the scattering matrix of the universal Hubbard model is 
performed and the BAE for some subsectors of the theory are determined. In section [6J a 
perturbative treatment a la Klein and Seitz [20] of the obtained Hubbard-like Hamiltonians 
is performed; second order and fourth order terms are presented. The last section is devoted 
to a short conclusion. Finally, we give in Appendix |A] some hints for progressing towards 
integrable bosonic Hubbard models, and expose in Appendix [B] a twisted version of XX and 
Hubbard models, leading to Hamiltonians that depend on phases that can be identified with 
a Aharonov-Bohm phase. 

2 Universal XX models 

We generalize the construction given in [16, 19,21] to the case of an arbitrary vector space 
V, possibly infinite dimensional. We will use the standard auxiliary space notation, i.e. 
to any operator A G End(V), we associate the operator Ai = A®I and A 2 = I <E> A in 
EndiV) <S> End(V). More generally, considering equalities in End(V)® k , Aj, j = 1, . . . , k, 
will act trivially in all spaces EndiV), but the j th one. 

To deal with superalgebras, we will also need a Z2 grading [.] on V, such that [v] = will 
be associated to bosonic states, v G Vo, and [v] = 1 to fermionic states, v G Vi. 

We will also assume the existence of a (super-)trace operator, defined on a subset of 
EndiV) and obeying cyclicity. When V is finite dimensional, dimV = n, EndiV) is just the 
algebra gl(n), so that the trace operator is the usual trace of n x n matrices. If V is graded 
and finite dimensional, one deals with the supertrace. When V is infinite dimensional, the 
definition of a trace operator is more delicate, and one needs to verify that it exists and is 
cyclic for the operators we use. We address this problem in appendix [A] 
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2.1 R-matrix 

To define the R-matrix of universal XX model, we need some preliminary notions. We define 
projectors 

tt: V -> W , 7f = I-n: W with V = W©W (2.1) 

and graded permutation operator 

lu<8>u -> (-l)NM v( g) U 
Note that in auxiliary space notation, the action of the (graded) permutation operator reads 

P\2U\v 2 = u 2 v 1 . (2.3) 
From these operators, one can construct an R-matrix 

R 12 (X) = S 12 P 12 + S 12 sinA + (I ® I - £i 2 ) Pia cos A (2.4) 
where £ i2 is built on the projection operators: 

£l2 = 7Tl 7f 2 + 7f 1 7T 2 (2.5) 



It is easy to show that £i 2 is also a projector in V <S> V: (£i2) 2 = £ 
Let us introduce the parity operator C: 

C = ix-tx. (2.6) 

It obeys C 2 = I and is related to the R-matrix through the equalities 

£12 = \{l - CiC 2 ) and I<g>I-E 12 = i(l + CiC 2 ) (2.7) 
that allow us to rewrite the .R-matrix as 
R{\) = cos(^) ( cos(^) P 12 + sin(~) I® l) - sin(~) C x C 2 ( sin(^) P 12 + cos(^) I ® l) . (2.8) 
One has 

Theorem 2.1 For a// spaces V and projectors n, the R-matrix \2.J$ satisfies the following 
properties: 

- Parity invariance: 

C 1 C 2 R 12 (\) = R 12 (X)C 1 C 2 (2.9) 

- Sign transformation: 

R 12 (-X) = C 1 R 12 (X)C 2 (2.10) 
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Symmetry: 
Unitarity: 
Regularity : 
Exchange relation: 



Pi 2 (A) = P 21 (A) (2.11) 
i?i 2 (A) R 2 i{-X) = (cos 2 A) 1(8) I (2.12) 
fl 12 (0) = P 12 (2.13) 



Rn(X) R2M = RM R21W (2.14) 
Yang-Baxter equation (YBE): 



H2J 



^12(^12) i?l3(Ai 3 ) #23(^23) — -^23(^23) -^13(^13) #12 (Ai 

where A^ = Aj — Aj. (2-15) 



- Decorated Yang-Baxter equation (dYBE): 

i?i 2 (A' 12 ) Ci -Ri3(Ai 3 ) -R23(A 23 ) = #23(^23) ^13(^13) Ci ^12(^12) 

Ay A, • A r (2.16) 

Proof: The proof is strictly similar to the one done in [16], the only needed relations being 

C 2 = I ; C\ Si 2 = Si 2 C\ = — Si 2 C2 = — C 2 S12 (2.17) 

and the relation (12. 7p . Let us also remark that this latter relation is equivalent to the 
relations 

2 Eij T, kj = + T, kj - E ik , V i, j, k with E ti = (2.18) 

without any reference to the projectors 7r and W. However, a detailed analysis of the relations 
(I2.18P shows that E i2 must be of the form (12.71) . up to conjugation. ■ 

Remark 2.1 When tt — or % — I, we get i?i 2 = cos(A) P12, which also obeys all the 
statements of theorem 12. 1\ but leads to trivial models. 

Lemma 2.2 If we denote by R^(X) the R-matrix built on it, we have 

R^(X) = Rfr\X). 

Proof: We have CM = -C^\ leading to the property E$ = sf 2 _7r) . ■ 
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2.2 Monodromy and transfer matrices 

From the R-matrix, one constructs the (L sites) monodromy matrix 

Aki...l>(A) = Roi(X) R 02 (X) ■ ■ ■ R 0L (X) (2.19) 

which obeys the relation 

Ro >{\ - fJ>) C 0< i... L >{X) Ay<i...i>(A*) = A'<i...L>(At) £o<i...l>(A) R 00 '(X - n) . (2.20) 

This relation allows us to construct an (L sites) integrable XX spin chain through the transfer 
matrix 

h... L {X) =tr £o<i... L >(A) =tr (i? 01 (A) i? 02 (A) • ■ ■ R 0L (X)) . (2.21) 

Indeed, when the trace operator is well-defined on the monodromy matrb|§, the relation 
f)2.20p implies that the transfer matrices for different values of the spectral parameter com- 
mute 

[*i..x(A),*i...l(m)] = 0. (2.22) 
Then, the XX-Hamiltonian is defined by 

^ = ti...L(0)" 1 ^(0). (2.23) 
Since the R-matrix is regular, H is local: 

L 

II X! ">■' 1 with H iJ+i = P M+^id+i ( 2 - 24 ) 
i=i 

where we have used periodic boundary conditions, i.e. identified the site L + l with the first 
site. 

2.3 Symmetry of universal XX models 

Proposition 2.3 Let us consider a universal XX model based on a vector space V, with 
projectors n : V — >• W and W : V — > W. For M G £W(W) © 2?nd(W), one has 

(Mi + M 2 ) i? i2 (A) = i? 12 (A) (Mi + M 2 ) . (2.25) 

As a consequence, the transfer matrix also has a symmetry (super) algebra End(W)®End(yV) , 
with generators given by 

M <L .. L> =Mi + M 2 + ... + M L , (2.26) 

where M e End(W) © End(W). The same is true for any Hamiltonian H built on the 
transfer matrix. 

2 For finite dimensional vector spaces, the trace operator is obviously always defined. For infinite dimen- 
sional spaces, one needs to be more careful: we will come back on this point in appendix 1X1 
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Proof: Starting from a general morphism M G End(V), a direct calculation shows that when 
M(W) C W and M(W) C W, we have Mtt = and MtT = f M so that fl^25|) holds. 
The above conditions are equivalent to M G i?nd(W) © £W(W). 

As far as the transfer matrix is concerned, the proof is the well-known, once (I2.25P holds. ■ 

Since the choice of the projector tx fixes W and W, the above procedure allows us to 
associate to any symmetry (super) algebra S a universal XX model possessing S as symmetry. 

The eigenstates of the transfer matrix will be also eigenstates of the Cartan genera- 
tors of the symmetry algebra. These generators are given by M aa , a = 1, . . . ,dimV = d 
(with possibly d = oo). The corresponding charges will be noted A = (Ai, . . . , A^). The 
charges (Ai, . . . , A r ), r = rank7r, correspond to End(W), while (A r+1 , . . . , A^) are associated 
to End(W). In the following, we will also need the fundamental weights 

A a = (0, 0,1,0,..., 0)*, a = l,...,dimV = d. (2.27) 

a-l 

3 Universal XX models based on gl(m\n) 
3.1 Hamiltonian and transfer matrix 

Projectors and i?-matrix: We apply the above construction to the case where V is the 
graded tensor product V = C m ' n , with possibly n = to encompass the case V = C m . In the 
following, we note s = n + m. 

The Z2 grading [.] is defined on indices j, such that [j] — 0, 1 < j < m, will be associated 
to bosons and [j] — 1, m+ 1 < j < tn + tl to fermions. Accordingly, the elementary matrices 
Eij (with 1 at position and elsewhere) will have grade [Eij] = [i] + [j]. 

To define the projectors ir and 7T, we introduce a subset 

j\fcz s = n z+ , 

and denote by A/" its complementary set, i.e. 

Af n AT = and A" U AT = Z s . 

We will also need the bosonic and fermionic 'components' of A/", 

Ao = A" n Z m Al = A" \ Ao with A/o U A/i = M . 

They are such that [j] = when j G Ao while [j] = 1 when j G A/i- 
To each set Af, one associates projectors 

7rW = J2 E jj > n = I s -7v = 7i iW) (3.1) 

Although these projectors depend on the set Af, we will drop the superscript (Af), keeping 
it only when several sets Af are considered. 
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From these projectors, one constructs the R-matrix according to the general formulas 
(12. 4p and (12. 5p . This .R-matrix obeys theorem I2.1[ with the parity matrix C : 

C = Y J E ]J -Y. E ^ = ^-*- ( 3 - 2 ) 



Monodromy matrix and Hamiltonian: From the R-matrix, one constructs the (L sites) 
monodromy and transfer matrices following the general procedure explained in section 12.21 
Then, the XX-Hamiltonian is defined by eq. ( 12.241) with two sites Hamiltonian 

Anticipating the Bethe ansatz analysis, one can see that this Hamiltonian describes, apart 
from the 'vacuum', m + n — 1 species of particles gathered into two subsets, so-called the 
'barred' a, b, . . . and 'unbarred' particles a, b, . . . corresponding to projectors W and n = I s —W 
respectively. The 'barred' particles move as hard-core particles while the 'unbarred' particles 
are displaced by the barred ones. This latter property is valid for a vacuum of 'unbarred' 
type: obviously, one has to reverse 'barred' and 'unbarred' particles if the vacuum is chosen 
of 'barred' type. 



Symmetry and number of models: Obviously, without any loss of generality, one can 
choose 

A" = {1,2, ...,r ; m + l,m + 2,...,m + ti} with r = \M \ , t x = |JVi| (3.4) 
A/" = {r + l,t + 2, . . . ,m; m + ti + l,m + ti + 2, . . . ,m + n = s} (3.5) 

From the property R^p(X) = R^\\) and the isomorphism gl(n\m) ~ gl(m\n), one can 
impose the inequalities 

. . . . ^ m + 1 . 
= I Ai | < — - — and m > n , 

leading to {[ m ^ ±1 ] + l) (min+1) different models, where we used the notation min = 
min(n, m) and max = max(n, m). 

The R-matrix admits a gl(m — to|n — ti) © gl(to\xi) symmetry superalgebra whose gen- 
erators have the form 

E jk ,j,keAf for gl(v Q \xi) , Q , 

Ejk , j,k£j\F for gl(m — x \n — Xi). 

As a consequence, the transfer matrix also admits gl(m — t |n — ti) © gl(x \xi) symmetry 
superalgebra, with generators given by 

M < i... £> =Mi + M 2 + ... + M L , (3.7) 

where M is one of the generators given in (13.61) . The same is true for any Hamiltonian H 
built on the transfer matrix. 
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Since the choice of the projector n fixes the values of to and ti, the above procedure 
allows us to associate to any symmetry (super) algebra S = gl(q\q') © gl(vx — q|tl — q') a 
generalized XX model possessing S as symmetry, provided the vector space V = C m ' n we 
start from is large enough (i.e. m > q and n > q' to get S). Conversely, from the vector 
space V = C m ' n , one can construct models possessing the symmetry: 

gl(q\c\') © gl(m— q|n— q') , q < m and q' < n (3.8) 



3.2 BAEs for universal XX models 

To get the BAEs of a model, one starts with a reference state, called the pseudo-vacuum, 
which is an eigenvector of the transfer matrix. The other states are constructed as 'pseudo- 
excitations' on this pseudo-vacuum. 



3.2.1 The pseudo- vacua sector 

The full space of states for the XX models is (V)® L : we consider here the subspace W vac = 
(VV)® L . In this subspace, the transfer matrix takes a simple form: 

txx(X) = (cosA) L P 1L P 2L . . .Pl-il + (sinA) L r with r = ranker (3.9) 

w„ ac 

One recognizes in txx(O) = exp(zp) the shift automorphism. The eigenvalues of p are the 
impulsions of the states. Note that the Hamiltonian Hxx = m (£)'(0) (given in ( 12.241) ) 
vanishes on this subspace. 

There are a priori r reference states 

Q a = (e a )' XlL , a = 1, . . . , r = rank7r (3.10) 

which have vanishing impulsion and charge LA a where A a is the fundamental weight given 
in (I2.27p . However, since the algebra S = End(W) © End(>V) is a symmetry of the model, 
one can restrict itself to highest weight vectors and get the remaining states through the 
action of the step generators of S. In fact, in (13.101) . there is a unique highest weight vector 

n 1 = ei®---®ei. (3.11) 

The other states in (13.101) can be obtained through iterative action of the symmetry gener- 
ators M Q) i: 

= (M a>1 ) L (3.12) 

In the following, we will take Qi is as the vacuum. The other states will be described as 
excitations above this vacuum, and we introduce for M indices bi, . . . , b M , and M positions 
Xi, . . . , xm the state: 

\{b} ; x >= eig - • - g ei ®e 6l © ei g ■ ■ ■ © ei ®e& 2 ©e x © • • •©ei©e bAr ©e 1 © • • •© e x (3.13) 

X\ — 1 £9— — 1 
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3.2.2 One excitation states 



We introduce 

L 



$i(p) = ^2e ipx \a,x> a = 2,...,rank(7r) =r (3.14) 

x=l 
L 

®l('P) = ^2e ipx \a,x> a = t+l,...,s = t + t = n + m (3.15) 



x=l 



where |a, x > is defined as in (13.131) . The indices a = 2, . . . , r correspond to the space W and 
the indices a — t+ 1, . . . ,S — t+t correspond to the space W. Through a direct calculation, 
it is easy to show that 

t(0)& a (p) = e ip ®l(p), u = a,a (3.16) 
H& a (p) = and H^l(p)=2cos(p)^l(p) (3.17) 



if p obeys the Bethe ansatz equation (BAE) 



(3.18) 



One can gather all these states into a single vector state. The set {a = 2, . . . , t; a = r + 
1, ...,£ = t + 1} is noted {j = 1, . . . , s — 1} where the first r — 1 indices are of type 'a' while 
the r last ones are of type 'a'. We introduce the elementary vectors Uj G C s_1 (with 1 at 
position j and elsewhere): they correspond to the 'small' chain of the nested Bethe ansatz. 
The vector state reads: 

s-l L s-1 

= j2 $j+i(p) Uj = J2 eipx \ x > with \ x >= Yl > u i ( 3 - 19 ) 

j=l x=l j=l 

Note that in \x >, \j, x > lies on the original 'big' chain (of length L), while Uj lies on a new 
'small' chain (here of length 1). The basic idea is to 'move' the action of the transfer matrix 
and symmetry generators from the 'big chain' to the 'small one'. Indeed, we have 

t(0)$ 1 (p) = e ip <$>\p) (3.20) 
H$\p) = D(p)^ 1 (p) with D(p) = 2cos(p)diag(0,...,0,l,...,l) (3.21) 

r-l x 

The matrix D(p) acts on the small chain (i.e. on the vectors Uj) while H was acting on the 
big chain (i.e. on the states \j, x >). In the same way, the charges of the states are given by 

M j+ld+1 <f>\p) = Eji&ip), j = l,...,5-l (3.22) 
Mn^fp) = (L-l)<$>\p) (3.23) 

where 6 End(C s_1 ), i,j > 1, (the elementary matrix with 1 at position and 
elsewhere) acts on the small chain. It corresponds to the generator of the symmetry 
generator My acting on the big chain. Mn (more precisely L — Mn) acts as a scalar and 
corresponds to the excitation number. 
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3.2.3 Two excitation states and scattering matrix 

We look for eigenstates $?-(pi,j? 2 ) describing two excitations of type % and j, and with 
impulsion p\ and p 2 . We gather these states into a single vector 

s-1 

® 2 (Pi,P2) = $ i+i,i+i(^i^2) Ui ® Uj - (3.24) 

defining a length 2 'small chain' (carried by the vectors Uj). The construction is done in the 
following way: 

$ 2 (pi,Pz) = {e ip ' x h-i®h-i + e i ^ x P 12 S 12 (p 1 ,p 2 )}\x 1 ,x 2 > (3.25) 

l<xt<X2<L 
s-1 

\xi,x 2 > = \i + + 1; Xx, x 2 > Ui<g> Uj (3.26) 

^i 2 (Pi,P2) = e~ ipi ^®7f+e iP2 7f®^-Pi2^®^ + 7f®7F) (3.27) 



where we have introduced p ■ x — p\X\ + p 2 x 2 and 7(p) • x = p 2 X\ +p±x 2 . P 12 is the (graded) 
permutation, 7f is the projector on W and n = I s _i — 7f. Since the scattering matrix S 
acts on the small chain, the projector n is the (lower rank) counter part in the small chain 
of the projector 7r (that acts in the big chain). We kept the same notation W for both of 
the projectors I s — n and I s _i — % because they are obviously isomorphic. This asymmetric 
situation is due to our choice of the vacuum, that belongs to W = 7r(V), not to W. 
The scattering matrix obeys Yang-Baxter equation and unitarity relation 

<Si 2 (pi,p 2 )S 13 (p u p 3 )S 23 (p 2 ,p 3 ) = S 23 (p 2 ,p 3 )S 13 ( Pu p 3 )S 12 ( Pu p 2 ), (3.28) 
S 12 (p 1 ,p 2 )S 21 (p 2 , Pl ) = I B . 1 0l s . ll (3.29) 
<S 2 i(pi,p 2 ) = Pi 2 S 12 (pi,p 2 ) P 12 , (3.30) 

while the braided S'-matrix <Si 2 (pi,p 2 ) = Pi 2 Si 2 (pi,p 2 ) (which appears in & 2 (pi,p 2 )) obeys 
braided Yang-Baxer equation and braided unitarity relation: 

S 23 (puP 2 )S 12 (p u p 3 )S 23 (p 2 ,p 3 ) = S 12 (p 2 ,p 3 )S 23 (p u p 3 )S 12 (p u p 2 ), (3.31) 
Si2{pi, P 2)S 12 {p 2 ,p 1 ) = I s _i®I s _i. (3.32) 

It is easy to show that 

t(0)<S> 2 ( Pl ,p 2 ) = e^ + ^$ 2 (p u p 2 ), (3.33) 
H$ 2 (p u p 2 ) = D(p 1 ,p 2 )^ 2 (p 1 ,p 2 ) with D(pi,p2) = D(pi)®I a _i+I,-i®D(p 2 ) 



if the BAEs 



e ip2L $ 2 (pi,p2) = S 12 ( Pl ,p 2 )® 2 (p u p 2 ) (3.34) 
^ L &(p u p 2 ) = S 21 (p 2 , Pl )$ 2 (p u p 2 ) (3.35) 
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are satisfied. We have used D{p) defined in (13.211) . leading to energies 0, 2cos(pi), 2cos(p 2 ) 
and 2cos(pi) + 2cos(p 2 ). The charges of the states are given by 



M j+1J+1 $ 2 ( Pl ,p 2 ) = (E jj ®I s _ 1 + I s _ 1 ®E jj }<f> 2 (p 1 ,p 2 ), j = l,...,s-l (3.36) 
MuK^bft) = (L-2)<$> 2 ( Pl ,p 2 ) (3.37) 

Again, the action of the Hamiltonian H and symmetry generators M.jj, j > 1 have been 
'moved' to matrices acting on the small chain. L — Mn is the excitation number. Since all 
the matrices are diagonal, we have indeed eigenvectors of the Hamiltonian and symmetry 
generators. 

Remark the property 

$ 2 (P2,Pi) = S 12 ( Pl ,p 2 y 1 $ 2 ( Pl ,p 2 ) = S 12 (p 2 , Pl ) $ 2 (pi,p 2 ) (3.38) 

that ensures that we can impose p\ < p 2 . 

The explicit form of the BAE depends on the type of excitation one considers. Looking 
at their projection on vectors Ui®Uj with i, j < r, one gets the BAE for type a, b excitations: 

e ipjL = u, j = 1,2 with u 2 = 1 (3.39) 

If one projects on ui (g) Uj with i,j > r, one gets the BAE for type a, b excitations: 

e ip > L = uj, 3 = 1,2 (3.40) 

If one projects on Ui (g) Uj with i < r and j > r, one gets the BAE for type a, a excitations: 

e i P i(L-i) = 1 and e i( Pl+P 2)L = 1 ( 3 _ 41 ) 

where p\ is attached to the type a excitation. 

If one projects on Ui®Uj with i > x and j < r, one gets the BAE for type a, a excitations: 

e i P 2(L-i) = 1 and e i( Pl+P 2)L = 1 ^ 342 ) 
where p 2 is attached to the type a excitation. 



3.2.4 M excitation states and BAEs 

We consider general states , $|jj(p), with M excitations of momenta p m , m = 1, . . . , M, 
gathered into a vector p. If M' is the total number of type 'unbarred' excitations, we note 
their corresponding momentum q n , n = 1, . . . , M', gathered in a vector q. In the same way, 
for M" = M — M' the total number of type 'barred' excitations, we noted q n , n = 1, . . . , M" 
their momentum, gathered in q. Hence we have 

{p!,p 2 , ...,p M } = { qi,Q2, ...,gM'}U{g 1 ,g 2) --- ) 9M"}- ( 3 - 43 ) 
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Then, the state $ A Jj(p) is characterized by: 



M excitations above the vacuum 

M 



M' 



M" 



n=l 



Momentum: \p\ = ^ p m = \q\ + \q\ = J^g n + g n 

m=l ra=l n 

M 

Charge w.r.t. the symmetry algebra: Q = Aj r 

M" 

Energy: E = ^ 2 cos (Qn) 



(3.44) 



m=l 



d = l 



where the weights Aj are given in (12.271) . 

All the states with M excitations can be gathered into a single vector 



*"(p) = £*&(p)«* 

{i} 



.72 



«,- M G C m ~ 1|n ® . . . ® C m ~ 1|n 



(3.45) 



describing a small chain of length M. Then, as for one and two excitation states, the action 
of the different integrals of motion can be 'transfered' from the original ('big') chain to the 
new 'small' chain: 



t(0) $ Af (p) 


= e i|p| $ M (p) 




(3.46) 


H$ M (p) 


= D(p)$ M (p) 

M 


M 


(3.47) 


Dip) 




_ x ®£>(p m ) <g> I s _i g> . . . <8> I s _! = V An(Pm) 


(3.48) 




i v 

m=1 m-l 


M-ra m=1 






= Ej_ij_i $ A/ (p) , j 

M 


= 2,...,r + f 

M 


(3.49) 


% 


= y i s _i ® . . . ® i a _ 


-i®Ejj ® I s -i ® . . . ® I s _i = ^ E^ 


(3.50) 




i' — 

m=1 m-l 


M-m m=1 




Mn $ A/ (p) 


= (L-M)$ A/ (p) . 




(3.51) 



In equalities ( 13 .461) - ( 13311) . the left-hand sides correspond to action on the original chain, 
while the right-hand sides corresponds to their counter-part on the 'small' chain. All the 
matrices in the r.h.s. are diagonal, and the projection of these r.h.s. on a generic state 
• • ® u jM reproduces the data (13.441) . 



Uj x <X> Uj 2 



The BAEs of the model take the form 

$ M (p) = S j+1J S j+2J . . . Suj S {j S 2j . . . Sj $ A/ (p) j = 1, . . . , M (3.52) 

where Sjk = Sjk(pj,Pk) is the two-body scattering matrix (13.271) acting in the spaces j and 
k of the tensor product explicited in (I3.45H . To compute them explicitly, we introduce the 
order -< defined by 

j + l^j + 2^...^M^1^2^...^j-l (3.53) 
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Any set of indices {ji,j2, ■ ■ ■ ,jn} ordered accordingly, ji -< j 2 -< . . . -< j n , will be noted 
\ji,32, ■ ■ ■ )in]-<- Then, from the form of Si 2 , one computes 

M 

Sj+2,3 ■ ■ ■ s 2j . . . Sj-ij — y~] (-i) n P331 ■ ■ ■ Pjjn { 

^ % n . . . n jn W kl ... 7f fcM _ 1 _ n exp (M - 1 - n) p^j 

M-l-n 

+ Tfj Tf h . . . Tf jn n kl ... h M -l-n eX P ( - * Pk ?) } ^ 3 - 54 ) 

1=1 

Above, the sum on {j}„ © {k} runs on partitions of [1,M] \ {j}, where the first set (of 
cardinality n) is ordered according to -<, {j}„ = [ji, j 2 , ■ ■ ■ ,j n }-<, while the second set {k} = 
{ki, k 2 , . . . , ku-i-n} is its complementary set. n k is the projector operator % in space k (of 
the 'small chain'). 

Applying (13. 54ft to $ A/ (p) leads to the BAEs for the XX model. To compute them, we 
remark that the operator Pjj 1 Pjj 2 . . . Pjj n corresponds to the cyclic permutation of the spaces 
Ji Ji) • • • i3n m the small chain. Its eigenvalues are (ui n ) k , k — 1, 2 . . . , n, where u> n = e 2j7r// ™. 
Moreover, this operator commutes with diagonal matrix within the brackets of eq. (13.541) 
and does not change the type of excitation. Thus, the BAEs take the form 

exp(zg n (L-M")) = (-l)*''" 1 (o^)" , n=l,2,...,M' (3.55) 
exp [iLq^j = {uj M "T exp , n = l,2,...,M" (3.56) 

M' 

with (ujm') M = 1 , (um") M = 1 and \q\ = q n 

n=l 

Remark that multiplying together all the BAEs one gets 

exp (iL \p\) = 1 . (3.57) 

4 Universal Hubbard models 

Starting with universal XX models, one can build universal Hubbard models, in the same way 
it has been done for usual and super Hubbard models [4, 16]. To simplify the presentation, 
we present the construction in the case of gl(m\n), but obviously the results are valid for any 
universal Hubbard model. 

4.1 R-matrices 

4.1.1 R-matrix for universal Hubbard models 

We start with the i?-matrices of two universal XX models, i?| 2 (A) and i?j 2 (A), leaving in two 
different sets of spaces that we label by f and j. Let us stress that the two XX models can 
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be based on two different (graded) vector spaces V] and Vj,, with two different projectors 7T| 
and 7T|, associated to the sets A/j and A/j_. 

The Hubbard model is constructed from the coupling of these two XX models. Its R- 
matrix reads: 

R&(\ lt A 2 ) = i4(A 12 ) 4,(A la ) + tanh C4) RUK2) C\ Ri 2 (X' l2 ) C{ (4.1) 

where again A i2 = Ai — A 2 and A' 12 = Ai + A 2 . The definition of the parameter h' 12 = 
h(Xi) + h(X 2 ) is given below. It is easy to show that this R-matrix is symmetric 

R\ l 2 (X 1 ,X 2 ) = R i 2 l(X 1 ,X 2 ), (4.2) 

regular 

/^(A 1 ,A 1 ) = Pg = P 1 t 2 j4 (4.3) 

and obeys the unitarity relation 

R\i(X 1 ,X 2 ) Rl{(X 2 , Ai) = ( cos 4 (A 12 ) - (^4v4 ^nh(ti 12 )) 2 ) l\ 2 s 4 

V. Vsm(A' 12 ) / / 

where I i2 = I <g) I (4.4) 



Property 4.1 When the function h(X) is given by 

sinh(2/i) = U sin(2A) (4.5) 



/or some (free) parameter U, the R-matrix ( 4-1 ) obeys YBE: 

^i2(Ai, A 2 ) #i 3 (Ai, A 3 ) i? 2 3(A 2 , A 3 ) = i?23(A2, A 3 ) i? 13 (Ai, A3) i? 12 (Ai, A 2 ) . (4.6) 
In that case, the coefficient in l[4-4\ ) can be rewritten as 

2 / 2 /tanh(/i 12 

where h± 2 = h(Xi) — h(X 2 ). 



^(X n)[ ^ { X n) - [ —^ )} (4.7) 



Proof: Again, as remarked in [16], the proof relies only on the properties (12.171) . (\2.7\i and 
follows the steps of the original proof by Shiroishi [22] , in the same way it has been done 
for algebras in [4]. Hence, the choice of the projector does not affect it. Moreover, it was 
already noticed in [4] that one can couple two XX models based on different gl(xn) algebras: 
this obviously extends to general (graded) vector spaces V. ■ 
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4.1.2 Gauge version of universal Hubbard models 

As for the usual Hubbard model, one can introduce a gauged version of the above R-matrix. 
It is defined by 

n 12 (x 1} x 2 ) = e ^^{ e iwj4 fl u (Ai)A2)e -iwicj e -ife4c» 

where hj = h(Xj) , j = l,2 (4.8) 

Being a gauged version of the previous R-matrix, 7^i2(Ai, A 2 ) also obeys YBE, is unitary and 
regular. This gauged version is used in usual Hubbard model to make contact between the 
above construction and the Hubbard R-matrix as it has been originally built by Shastry. 

4.2 Monodromy matrices, transfer matrices and Hamiltonians 

We remind that for given vector spaces Vj. and Vf, the different possible projectors 7Tj_ and 
7T| give different R-matrices with, as we shall see, a different symmetry (super) algebra. 
We consider the 'reduced' monodromy matrix 

(A) = TZ abl (A, 0) . . . TZ abL (A, 0) (4.9) 

and, when the trace is well-defined, its transfer matrix 

t(X) = tr a L a<bl ... bL> (X) 

Then, one gets 

[H,t(X)} = 0, VA , for H = H(0) =*(0)~V(0) (4.10) 

This 'reduced' monodromy matrix is just the one used to define the Hubbard model; one 
can compute 

^ 12 ^'°) = ^7jTli l (h)Rl 2 (X)R{ 2 (X)li l (h) (4.11) 
cosn^n.) 

where 

I^(h) = cosh(^) I <g> I + sinh(^) C\ C{ (4.12) 
The explicit form of the Hubbard Hamiltonian reads 

L 

H = Y,H hJ+1 (4.13) 
j"=i 

with 

",j • i P] d+1 + Sj ii+1 P} J+1 + UC]C} (4.14) 

where we have used periodic boundary conditions. One can see that the kinetics is dictated 
by the XX models: barred particles moves 'almost freely' with the noticeable exception 
that a} and F (or and Ir) cannot cross. Unbarred particles of type up (resp. down) 
are displaced by barred particles of same type. There is interaction only between 'up' and 
'down' particles, and the sign of the interaction depends on their 'bar' or 'unbar' type. 
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4.3 Symmetries 

We generalize the results obtained for su(m) Hubbard models (see for instance [4,11]) and 
gl(m\n) Hubbard models [16]. 

Proposition 4.2 The transfer matrix of generalized Hubbard models admits as symmetry 
(super) algebra 

EndiVl) © Endiyl) © End(V^) © End(y[) , 

each of the EndiV^) © EndiVl), e =f, I corresponding to the symmetry of one XX model. 
As a consequence this symmetry is also valid for the Hubbard Hamiltonian. 

Proof: To prove this symmetry, it is sufficient to remark that 

MC = CM (4.15) 

where M = M T + M l and M e G End(V^) © End(Vl), e =], |. Thus, one gets 

[R 12 {\, 0) , Mj + M|] = = [R 12 {\, 0) , M| + Mi] (4.16) 

where i?i2(A, 0) is the f?-matrix of the universal Hubbard model. 

As far as Hamiltonians and transfer matrices are concerned, the generators of the sym- 
metry have the form 

L L 

M T = J^M] and M j = ^M] (4.17) 



The eigenstates of the transfer matrix will be also eigenstate of the Cartan generators M^ Q 
a = 1, . . . , dimV 6 = d e , e =f , |. The corresponding charges will be noted A e = (\\ , . . . , X £ d ) 

5 BAE for universal Hubbard models 

We follow the same steps as in section 13.21 

5.1 Scattering matrix 
5.1.1 Pseudo- vacua sector 

The full space of states is now (V T © V J -) (g)L , and we consider the subspace W mc = (W T £> 
W^)® L . In this subspace, the Hubbard transfer matrix takes a factorized form: 

t(\) =(l + tanh(/ i )) L 4 x (A) t xx (X) (5.1 

Vvva 



Wiac 



where h(X) has been given in (14. 5p . The eigenstates of this sector also take a factorized form 

^fmM) = (5.2) 
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with eigenvalues 

£(p,p';A) = (l + tanh(/i)) L ((cosA) L e i|p| + (sinA) L t T ) ((cos A) L e* |p '' + (sin A) L . (5.4) 
Above, we have introduced v = rank(?f T ) and tr = rank^). The charges of the states read 



771=1 



M Af' 



M? *OhS'}(P> P') = A? *OhS'}(P» PO with At = (L - M) A! + A L (5-5) 

m=l 
M' 

^un(P>p') = Al ^un(P>p') with A^ = (L-JiOA} + X)Ai. (5-6) 

Their momentum is given by 

P^-UfMP') = nn^(p,p';0)$;-{;, } (p,p') = {J2p m + Y,p m )^uf^p') 

771=1 771=1 

= (IpI + Ip'I)^}^^') ( 5 - 7 ) 

and their energy reads 

Tj-foMM' ( n n i\ _ [d\£(PiP j ^)]a=0 frMM' / /% _ rr r ^Af Af' /_ _/\ /c o\ 

^^oMn^'P) - — s(p,p>; o) — - UL ®{j},myP>P) y 5 - 8 ) 

5.1.2 General excitations 

Now, we perform general excitations above the vacuum fij ® We note 

s e = n £ + m £ = r £ + r e , r e = rank(7r £ ) , F = rank(7f £ ) , e =T, J. (5-9) 

3=5^+S l . (5.10) 

We will have four types of excitations: 

'''unbarred' of type up: (a, |) = j + 1 , 1 < j < r T — 1 

'barred' of type up: (a, |) = j + 1 , r T < j < s T - 1 1 

'unbarred' of type down: (a, j) = j + 1 — s T , s T < j < s T + — 2 

k 'barred' of type down: (a, j) = j + 1 — s T , s T + — 1 < j < s — 2 

The set {(a, e)} corresponds to the space W e , e =f, I- The set {(a,e)} corresponds to the 
space W e without the index 1 (that is associated to the vacuum): in the following, we will 

o 

note this reduced space W £ . 
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One excitation For the states with one excitation, one has just to mimick what has been 
done in section 13.2.21 

L 

®l,e(p) — ^2 e%px l a > £ : x > a = 2, . . . , rank(7r e ) = x £ 

x f , e=U (5-12) 

®1, s (p) = ^2 e%PX I"' £ ' x > a = x £ + I, . . . ,s £ = x £ + x £ 

x=l 

Through a direct calculation, it is easy to show that 
((0)< £ (p)=e !p < £ (?), a = a-a 

H<$>\ E {p) = UL and H ^(p) = (2 cos(p) + U(L - 2)) ' 

if p obeys the Bethe ansatz equation (BAE) 

e ipL = 1 (5.14) 

Again, one can gather all these states into a single vector state. The labelling of excitations 
is done as explained in (15.111) : 

{a = 2,...,r T ; a = t T + 1, . . . ,s T = r T +r T } -> {j = 1, . . . , s T - 1} (5.15) 
{a = 2,...,i 4 ; = ^ + 1,...,^ = ^+^} -> {j = s T ,...,5-2} (5.16) 

where the first — 1 indices are of type 'a, f, the next ft are of type 'a, f, and so one. 
We introduce the elementary vectors Uj G C s ~ 2 (with 1 at position j and elsewhere) 
corresponding to the 'small' chain of the nested Bethe ansatz. The vector state reads: 

8-2 0-2 

Note that in |x >, \j,x > lies on the original 'big' chain (of length L), while Uj lies on a 
new 'small' chain (here of length 1). As in section [3.2.21 we 'move' the action of the transfer 
matrix and symmetry generators from the 'big chain' to the 'small one'. We get 

t(0)$ 1 (p) = e ip <$>\p) (5.18) 
H<$>\p) = V{p)&(p) = ((2cos(p) -2U)D + E/LI 8 _ 2 ) <$>\p) (5.19) 

D = diag(o,...,0,l,...,l,0,...,0,l,...,l) (5.20) 




t 



T-l 



The matrix T>(p) acts on the small chain (i.e. on the vectors u 3 ) while H was acting on the 
big chain (i.e. on the states \ j, x >). In the same way, the charges of the states are given by 

M j+1J+1 &(p) = Eji&ip), j = 1,... jJ6 -2 (5.21) 
<5>\p) = (L-D e )^(p), e=U (5.22) 

D T = diag(l,...,l,0,...,0) and D i = diag(o, . . . ,0, 1, . . . , l) (5.23) 



In 



sT-l si—1 sT-1 0-L-l 
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where Eij G End(C s ~ 2 ), i,j > 1, (the elementary matrix with 1 at position and 
elsewhere) acts on the small chain. It corresponds to the generator of the symmetry 
generator acting on the big chain. (L — M. £ u ) corresponds to the excitation number for 
e particles (e =], |). 



Two excitations For more than one excitation, a new effect appears with respect to the 
XX models: there can be two excitations at the same site (provided there are of j and | 
type). To take it into account, we perform a change of basis on the states and define: 

S-2 s-2 

x 2 >= ^2 I* + 1>.7 + Xl > X2 > Ui ® u 3 (5-24) 



t=i j=i 



with the convention that 



\i + 1, j + 1; x, x >= < 



(0 i,j<s T -l 

| \i + l,x > ®\j + l,x > z < s T - 1 < j 

+ 1,£ > l,x > j < s T - 1 < z 





s T - 1 < i, j 



Then, the eigenstates are gathered into a vector 



s-2 



i=i j=i 

We have 

$ 2 (Pi,P 2 ) = ]T {e ip -* I s _2® 1,-2 + e i ' r(p) - a: Pi2«S 12 (p 1 ,p 2 ) 

l<a:i<X2<-I/ 

The scattering matrix is given by 



\Xi, x 2 > 



S? 2 E {P1,P2) 



e~ vpi ¥ ® W £ + e ip2 lf £ 



^ <8> (rr l + 7T 1 ) + (n l + Tf 1 ) <g> £ T + 7T l (g> 7f T + 7f T <g> ^ 
(t( Pi ,P 2 ) I s _ 2 ® I s _ 2 + R(pi,p 2 ) Pis) ® Tf 1 + H 4 ® 7f T 

sin(pi) - sin(p 2 ) 
sin(pi) — sin(j9 2 ) — 2iU 
2iU 



S 



Sl 2 {pi,P2) 

T(pi,p 2 ) 



sin(pi) — sin(p 2 ) — 2iU 



T(pi,p 2 ) ~ 1 



(5.25) 



(5.26) 
(5.27) 

(5.28) 

(5.29) 

(5.30) 
(5.31) 

(5.32) 
(5.33) 



where rc £ (resp. 7f £ ) is the projector on W £ (resp. W £ ), e =f, J,. 

One recognizes in 5^ e (pi,p 2 ) the scattering matrix of an XX model in the l e subsector' 
(e =t)l)- They correspond to the only part of S which acts non trivially in the jf and 
|| sectors. The remaining part (acting in the || and || sectors) have been divided into a 
part acting only in the 'bar sector' (the S H matrix, of Heisenberg type) and the rest (the 
matrix) . 
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Comparison with usual Hubbard model: the parts Si 2 e {pi,P2) and S^ipiiP?) in the 
scattering matrix are just generalizations of the Hubbard scattering matrix to higher dimen- 
sional case. Note however the projections appearing in these scattering matrices, that are 
new w.r.t. the usual Hubbard model: we will comment on this point in section 15.2.41 The 
part S* is completely new: it introduces new physical effects that were not seen in Hubbard, 
due to the 'small size' of its vector space. 

5.2 BAEs: a first account 

Once the scattering matrix of the universal Hubbard model is known, the technique to 
obtain the transfer matrix eigenvalues and the BAEs is a priori known, see section 15.2.11 
below. However, if the eigenvalues are easy to deduce, the determination of the precise form 
of the BAEs is a more delicate problem. Here, we compute them for some subsectors of the 
theory, leaving the determination of their complete form for a further publication. 

5.2.1 M excitations 

We consider a general state with M excitations, that divides into excitations of type a 
in the "f sector', excitations of type a in the "f sector', excitations of type a in the 
'J, sector', and excitations of type a in the 'j. sector'. The construction follows the line 
of section 13.2.41 with the noticeable exception that there can be j and I excitations at the 
same site. To take this fact into account, we introduce: 

| j m \ Xm > if all XmS are different 

(5.34) 

if at least three x m } s are equal 

'0 iij,j'<s^-l 

if j < s T - 1 < / 

(5.35) 

if f < s T - 1 < j 

if s T - 1 <j,f 

Then, the BAEs take the form 

$ M (p) = S j+1J S j+2J . . . S Mj Sy S 2J . . . S s ,. ; $ A/ (p) j = 1, . . . , M (5.36) 

In the following, we examine the BAEs in subsectors that are related to different types 
of excitations: the two XX-type subsectors, where excitations are only of type j or only of 
type | ; the 'unbarred subsector', where excitations are only of type 'unbarred' (f or j), and 
the Hubbard-type subsector, where excitations are only of type 'bar' (f or |). 

5.2.2 BAEs for the XX-type subsector 

We introduce the projectors on the ' j sector' and 'I sector' 

IT = n £ +W £ , e=T,l . (5.37) 



\{ihx> 



\ 


9) M 






\ 


>?) M 




&) m=l 
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It is easy to see that 

it ® u e s 12 ( Pl , P2 ) = sX £ ( Pl , P2 ) it ® rr e =t, | (5.38) 

so that multiplying the BAEs from the left, by (II £ )® M , one recovers the BAEs of the XX 
models: 

e - ^j+2,j ■ ■ ■ &Mj *>\i ■ ■ ■ *t yp) 

j = l,...,M ; M l = M l = (5.39) 

e q> l (p) - & j+1J & j+2 j . . . & Mj by j> 2j . . . . j g> T (p) 

j = l,...,M ; M T = M T = (5.40) 

$fb) = ( n £ ® it ® . . . n £ ) $ M ( P ) , £ =T,| (5.4i) 

M 

These BAEs corresponds to subsectors where excitations of only f or only J, types are con- 
sidered. They are of the same form that the XX models: 

exp (i g n (L - M")) = (-l)^'- 1 (w M ,) n , n = 1, M' with (o; M M ' = 1 (5-42) 

exp (iLq n ] = (uo M ») n exp ( - z |qrA , n = 1, M" with {u M ») M " = 1 ■ (5.43) 

5.2.3 BAEs for the 'unbarred sector' 

We consider the 'unbarred subsector', i.e. states with unbarred excitations (of type j or |) 
only. The corresponding projector is 



n = £ T + ^ ; n 12 = n®n ; n L „ M = n®n®. . .®n®n (5.44) 

From the property 

n 1 ... M S 12 (p 1 ,p 2 ) = Sg(p 1 ,p 2 )IL 1 „ M (5.45) 

e=T.J. 

it is easy to see that the calculation is very similar to the XX case, with 7^, 7^ playing the 
role of 7T, 7f of section 13.2.41 Using the same notations, one gets 

M 

Qun cun Qun cun cun cun \ ^ \ ^ / -i \n p TD TD J 

e Y I,; i+2j • ' ' ^Mj «>y ^2/ • • • 2^ 2^ ' U ' -//':• • • • 1 \ 

-1 — 71 J 

Since in this sector the BAEs take the form 

6 ^un\P) — °j+2,j ■ ■ ■ <~>Mj °lj °2j ■ ■ -dj-lj ^un\P) 

C(p) = n L .. M $ M ( P ) 



n=0 U\l 



{jK®{k} 



Ot Of Of 01 01 OIOI 1 Of ot 

7T 7T • . . . TC ■ 7Tt - - - 7Tt f , + 7T| 7H. . . . 7T} 7T.L . . . 7rl 

J Jl Jn Kl Km— 1 — n J Jl J" n-1 n>M- 
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we obtain 

exp(iq n L) = (-l)^- 1 (o; MT ) n , n = 1, 2, . . . , M T with (w mT ) mT = 1, (5.47) 
exp^i) = (^Mi) n , ra = 1, 2, . . . , JW- 1 with (w Mi ) Mi = 1 (5.48) 

where is the number of | excitations, and = M — is the number of j excitations. 
We noted q n , n — 1, 2, . . . , M^, the momenta of the j excitations and g^, n = 1, 2, . . . , M^-, 
the momenta of the j excitations. 

Remark that the two series of BAEs (I5.47P and (15.481) are decoupled, and correspond to 
the 'unbarred sector' of each of the underlying XX models. They can be obtained separately 
using the projectors 

n i...M = ^ T ® ^ ® ■ ■ ■ ® ^ ® k\ or U[ M = ® tc 1 <g> . . . ® (8) (5.49) 

Af M 

but the present calculation shows that they are complete in this subsector. 
5.2.4 The 'bar subsector' 

Following the same lines as in the previous sections, one can consider the 'bar subsector', 
i.e. states with a<\ and a± excitations only. The corresponding projector is 

n = w T + ¥ and n L .. M = n ® n ® . . . ® n . (5.50) 

M 

In that case, one has 

ni...M«5i 2 (pi,p 2 ) = s 12 {pi,P2) ni...M = s 12 {pupi) n L .. M 

S 12 {puP 2 ) = fT(p 1 ,p 2 )I s _ 2 ®I s _ 2 + J R(p 1 ,p 2 )P 12 N ) f 7f T ® Tf- 1 + VfJ- ® 7f T ) 



-F i2 (^7f T ® 7f T + ® TfJ-J (5.51) 

One could be tempted to reccognize in <S 12 , the scattering matrix of a generalized XXX 
model. Indeed, specifying the f or j type only (whatever the indices a,b,c, . . . are), one gets 
on a state with two excitations: 

s_i 2 ( Pl ,P2)m'> = -m> (5.52) 
<s 12 ( Pl , P2 )iu> = T( Pl , P2 )\n>+R(pi,P2)m> (5.53) 

S 12 ( Pl ,p 2 )\U> = R(px,p 2 )\n>+T(pi,p 2 )m> (5.54) 
«5i 2 (pi,p 2 ) HI' > = -||'|> (5.55) 



In the case of Hubbard model, where there is only one type of f excitation and one type of 
j. excitation, one has exactly the scattering matrix of the XXX models. This allowed the 
calculation of BAEs of the Hubbard model. However , if there is more than one type of j or 
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I excitation, this is not the case anymore. For instance, for two types of j excitations (say 
a and b), eq. (I5.52p corresponds to 



5 12 (pi,jP2) \a T;o T' > = -|at'; a T> (5-56) 

Si 2 (pi,p 2 )|aT;&T'> = -|&T';aT> (5.57) 

S 12 {pi,P2)\bV,aV > = -\aY;b^> (5.58) 

S 12 (p llP2 )\bV,bY> = -|&T';&T> (5.59) 

while a generalized XXX model would act as 

Sizipu&z) \a T;a f > = -|at';at> (5.60) 

«5i 2 (pi,p 2 ) |a t;& t > = r (Pi>P2) |a t;& t > +R(pi,P2) \b t;a t > (5-61) 

5 12 (p 1 ,p 2 )|6T;at> = T(p 1 ,p 2 )\bV,a J [>+R(pi,p 2 )\aV,b J [> (5.62) 

5i 2 (pi,p 2 )|6T;6T'> = -|&T';&T> (5.63) 

This difference just prevents to perform a nesting in the same way it is done for Hubbard. 
Note that <Si 2 still obeys 

S 12 (p,p) = -Pv (5.64) 



so that one can define an integrable spin chain associated to the nesting in the usual way. 
However, the exact form of the BAEs for this new chain is not known yet. The same is true 
for the general BAEs of the universal Hubbard model. We will come back on this point in 
a further work [23]. 

6 Perturbative expansion of the Hubbard-like Hamil- 
tonian 

We expand the Hamiltonian ( 14.131) and (14.141) in the inverse coupling jj. That expansion 
has been used in [13] to match the SU(2) dilatation operator with the effective Hamiltonian 
of the Hubbard model. The system was taken at half-filling to guarantee the required spin 
chain behaviour. 

Being ultralocal, the potential term U^2j CjCj is separately diagonalizable on each site 
with eigenvalues rkU. Indeed, they are obtained from the property C 2 = I ( 12. 6p . 

1* I 

The ground state has eigenvalue —LU and can be obtained if the condition Cj Cj = — 1 is 
realised on each site. This is equivalent to demand eigenvalue 1 on each site for the (one-site) 
projector 

1-6^ = JizMi j = %] + ^ _ 2 , H = _ ^ (6 . 1( 
or the global projector 

no=n«-*i) 2 =nw-^) 2 = n 5- < 6 - 2 ) 

3 j 
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We observe that it projects on the subspace where, on each site, one and only one projector 
among 7fj , 7fj has nonzero action. This means that only the following subspaces survive 

n : or W ®W [ (6.3) 

3 3 

Namely we demand that on each site there is one barred particle: double or empty occupan- 
cies of barred particles are prohibited. This may be possible only if the system has precisely 
L barred particles out of the 2L permitted ones. We say that the system is half-filled and 
we will assume this condition to perform the perturbative calculations. 

It is useful to compare with the ordinary Hubbard model, where the algebra is realised 
in terms of fermionic oscillators c a j, c^. ■, satisfying {c a j, c£ ■} = 1 (a =f , J,). There, 
the projector 7fJ is equal to the number operator n a j = c a jC aj j, so in the present general 
formalism, a vector of W a corresponds to an electron and a vector of W CT corresponds to a 
vacancy (in the Hubbard model we have tt J = 1 — n a j). 

We follow the method introduced by Klein and Seitz [20]. With reference to the Hamil- 
tonian (14.141) and in complete analogy with previous cases [16] , we define a hopping operator 
by 

Intuitively we associate its action to the move of a barred particle from site j to site i. We 
define hermitian conjugation as super-transposition (because our operators are real) so that 
we have 

4 = x ^ ( 6 - 4 ) 

The two-sites Hamiltonian (14. 14ft takes the form 

X jJ+1 + X j+1J + UC]q (6.5) 
and is obviously self-adjoint. The perturbing term is 

r = E(x w+1 + x i+1J ). (6.6) 

3 

The action of Xij on the vector spaces defined in (12.11) is easily described by observing the 
projectors in the second line of (j6.4p . We initially focus on the effect on site i: 



1) 


W T 


® w l 
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> 






2) 


W T 
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i 




3) 


W T 


®w l 
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Xij 
> 
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4) 


W T 


<8> W l 


Xij 
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W T 





(6.7) 
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On site j the description is the complementary one, namely: 

1) W^(g)W l 
j 



J (6.8) 

3) W T ®W — W T ®W A 



2) W'®W l — ^ W T ®W* © W 1 0W l 
j 

4) W T ®W j W T ®W j 

It is clear that the domain and the codomain of X^ are always disjoint and only a "two-fold" 
action can make them the same. We state this as a theorem. 

Theorem 6.1 Given a site i and an initial configuration choosen among the four listed in 
(6.1), the product of an odd number of operators X^ or Xji acting on i with possibly different 
j cannot return to the same initial configuration. 

Proof: The proof is a trivial application of the rules in (16.71 and 16.81) ■ 



A number of corollaries follow from it and from the action of the projector Ho: 

Xfj n = (6.9) 

(i-n ) x i:j x^ n = o 
x j-i,j x j+ij n = o 

n T n II = if n = odd and L > n 

The condition L > n in (16. 9p is extremely important. If it is removed, new terms known 
as demi-wrapping L = n and wrapping L < n behaviours occur because the perturbative 
interaction (16.61) circulates all around the periodic chain. These terms will not be evaluated 
here. 

According to Klein and Seitz [20] , the effective Hamiltonian for the states originated from 
Ho is 

H«=±H$ + ±H<$ + ... (6.10) 



where 



s = (i-u )—l— (i-no) (6.ii) 
= n TSTn = ~n t 2 n (6.12) 
= ^n T 2 ST 2 n + ^n T 2 n T 2 no (6.13) 

These expressions can be worked out following Klein and Seitz [20], with a long but simple 
calculation that is not shown here, using the properties (16. 9ft . 
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We point out another property that allows to simplify the expressions: the projector Ilo 
makes redundant one among the barred and the non-barred projectors. Indeed, we can write 
the identity on site j as 7rJ + 7fJ = 1, therefore 

7rjn = ttJ {ttJ° + 7fJ a )U = TvJnj^Uo = WJ a U (6. 14) 

and we can write all the expressions using only the non-barred projectors. 

6.1 Second order Hamiltonian 

A direct calculation shows that for L > 2 the second order effective Hamiltonian is 



The structure of the two-sites Hamiltonian P cff 1 2 can be described in the following way. 
The projector n allows states of the form (16.31) . so we start observing that the Hamiltonian 
vanishes on states of the following form 

W^®W l ® W^W 1 or W ®W [ ® W®W l (6.16) 

12 12 

because the projectors in (16.151) require orthogonal subspaces on different sites for the same 
type (e.g. up) of vectors. For example, this means that a state ®W <g> ® tr is killed 

1 2 

by the two-sites Hamiltonian. We are left with states of the form 

W^W 1 ® W T ®W l or W ®W l ® W^®W l (6.17) 

12 12 

on which the parenthesis of projectors [irj 7ij +1 + 7rj 7r| +1 ) acts as the identity. A state in 
one of the spaces (16.171) is respectively of the form 

v®W <g> v®w , v0w <g> v®w (6.18) 

12 12 

on which we have respectively 

(l + Pj,, tP/,, ,) u(8)uJ ® U®iy = UOw + (-1)(M+N)([«]+M) ® 

V JJ-t-J- JJ+J-/ 12 12 12 

(6.19) 

(l + PLiPL,) v®w v®w = v®w u0w + (-l)(H+[»])([^+M) ^ 

\ J,J-r j,jt i 1 2 1 2 12 

In matricial form, the two-sites Hamiltonian has one of the two block-diagonal structures 

B - = (-1 ~l) ° r B+ ={\ l) ' (6 ' 20) 
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all other entries being zero [16]. Both the blocks have eigenvalues and 2. The multiplic- 
ity depends on the actual model under examination. In the Hubbard model, the effective 
Hamiltonian acting on the singly occupied states reduces to the block B_ only 



(2) 
eff 



1-4S 



/0 0\ 

1-10 
0-110 

\0 0/ 



(6.21) 



where Sj = (S x , S y , S z ) on site j are the spin vectors of the Heisenberg model. 

6.2 Fourth order Hamiltonian 

The fourth order Hamiltonian 



(4) 
off 



£*£W a = (6-22) 

3 

3^ £ + 2 - P 2i+l- P ij+l + ^i+lJ+2^/+lJ+2^ii+l jP ii+l) ( 7r j 7r ]+l 7r j+2 + 7r ] 7r j+l 7r j+2) + 

+ (l + 2P J T +1J+2 Pj +lj - +2 + P iJ+1 P ij+1 P i+lij+2 -P, +lij+2 ) ( 7r ) 7r )+l 7r ]+2 + 7r j 7r "j+l 7r "j+2) + 



-2(2 + / 'j.j . I / % . I + P/+1J+2-P/+1J+2) ( 7r j 7r j+l 7I "j+2 + 7r ) 7r j+l 7r j+2) 



n, 



is composed by a three-sites Hamiltonian density. It acts generically on states 7r J 1 7r^+i 7r J+2 ) 
where e = ±1 indicates respectively | and [. It cannot mix states with different values of 
£i + £ 2 + £3; if that sum is ±3, its action is zero. The other two possible values are ±1, on 
which it acts independently. 

The second order Hamiltonian can be put in a three-sites density form by averaging on 
neighboring sites 



(2) 



H. 



(2) 



2 v eff ~ J1 eff 3+1,3+1) 

therefore we can evaluate the eigenvalues of a three-sites Hamiltonian formed by 



K 



eff ij+lj+2 



-(h (2) +h (2) 



(4) 



^ jj2 MA eS j,j+l,j+2 



(6.23) 



(6.24) 



by action on states as in (I6.19p . The possible eigenvalues are (up to corrections of order 
U~ 4 ) the same in all cases 



Eigen(ft cff 



3,3+1,3+2) 



0,1,3 1 + 



IQU 2 



(6.25) 



with multiplicities that depend on the specific model under consideration. 
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7 Conclusion 



We have defined in a very general way Hubbard models with arbitrary symmetry. The basic 
ingredients are a vector space, which defines the representation space on each site, and two 
projectors, which separate the particles into two classes that behave in a different way. The 
scattering matrix, as well as the energies, have been computed on very general ground. The 
general form of Bethe ansatz equations remains to be computed, although some of them 
are given here. Of course applications of these models to condensed matter physics and/or 
AdS / CFT correspondence are of first importance. In this regard, we have given in appendices 
a general procedure to include a Aharonov-Bohm phase and some hints towards the definition 
of integrable bosonic Hubbard models. This latter feature is especially crucial, in particular 
for applications in string theory regarding the AdS/CFT correspondence. Indeed, dealing 
with models with psw(4|4) symmetry, one is lead to consider some subsectors of the theory 
containing both fermionic and bosonic particles. The study of condensed matter models 
with bosonic content on an integrable point of view also requires progress in this direction. 
Finally, let us stress that the input of boundaries, which play a great role in condensed matter 
models, may be a worthwhile extension of this work, see e.g. [24] for the one-dimensional 
Hubbard model with integrable boundaries. 



A Towards integrable bosonic Hubbard models 

Our construction can in principle be applied to infinite dimensional vector spaces, leading 
to possible bosonic integrable Hubbard models. However, for such a purpose, one needs to 
construct a trace operator on the space V, a task that is not guaranteed when V is infinite 
dimensional. To simplify the presentation, we work on XX models, but the procedure leading 
to Hubbard models can be applied in the same way we did for finite-dimensional vector 
spaces. 



A.l R- matrices associated to Fock space 

To illustrate the problems encountered in the infite dimensional case, we focus on the case 

where V is the Fock space based on oscillators (b,b>) with [&, ^] = 1. We denote by 

|0 > the Fock vacuum, and \n >= (b^) n |0 >, n — 0, 1, 2, . . ., and by N = b the number 
operator: 

N \n >= n\n> . (A.l) 

The dual vectors are noted < n\: 

< m\ . \n >=< m\n >= 5 n , m . (A. 2) 

We also introduce the subsets 

F N ={\n>, n = 0,1,2, ...,n\ C F, N = 0,1,2,... (A.3) 
The permutation operator is given by 

oo 

P12 = ^2 \n >< m\ <g> \m >< n\ (A.4) 

n,m=0 



28 



We present two examples of projectors. 



Even-odd projectors: One can choose as projectors 

n ev = i(l + (-1)") and T ev = I — n ev = n odd = i(l - (-1)*) , (A.5) 

which project on even and odd particle number eigenspaces T ev and T dd- They obviously 
commute with N. The parity operator is C = (— 1)^. Then, the corresponding XX i?-matrix 
reads 

R(X) = cos(^) { cos(^) P 12 + sin(~) I® ij - sin(~) (-1)^+^ { sin(^) P 12 + cos(^) I ® l} 
where N ± = N ® I and N 2 = I ® N. 

Small modes projector: For any number I G Z + , one can also take as projectors 

i 

7r- e = ^2 \ n >< n \ and lr- e = I-n- e = tv >£ = ^2\n><n\. (A. 7) 

n=0 n>£ 

They also commute with N. The parity operator reads 

& oo 



c = \ n >< n \ - i n >< n i • ( A - 8 ) 



n=0 n=l+l 

Properties of the .R-matrix: The different types of projectors lead to different types of 
i?-matrices. From the general treatment done in section 12. f I one already knows that these 
i?-matrices obey the theorem 12 .11 in particular the Yang-Baxter equation. Then, one directly 
constructs a monodromy matrix 

£o<i...l>(A) = j Roi(A) j Ro 2 (A)--- j R l(A) (A.9) 

which obeys the relation 

^oo'(A - £t) Aki...l>(A) £ '<i...i>(A t ) = £o'<i...£>(/-0 £o<i...l>(A) Rq >(X - fi) . (A.10) 

However, to get a transfer matrix leading to an integrable model, one first needs to define 
the trace operator: we discuss it in the next section. 

We recall that the same construction is valid for the Hubbard f?-matrices that one could 
build by coupling two XX models, as done in section 14. 1. 11 

A. 2 Trace operator and transfer matrix 

One can define a trace operator for operators O such that [O , N] = in the following way. 
Since such operators preserve the subsets J 7 ^, which are finite dimensional spaces, their trace 
on jF/v 

N 



tr N = Y^ < n\0\n> (A.ll] 



n=0 
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is well-defined and cyclical The trace on T is then defined by the inductive limit 



trO = lim —tr N 0. (A. 12) 



It is cyclic and such that 



M=l ; trir ev = trn ev = - ; trvr- £ = and trvf^=l. (A.13) 

In the same way, the permutation operator commutes with Ni + N 2 , so that £r 12 Pi2 is well- 
defined and cyclic. However, to define the transfer matrix, we need to use the partial trace 
tr±Pi2 which is ill-defined because P12 does not commute with Ni. For instance, it is easy 
to see 

tr Nl (P12 P13) ± tr Nl (P 13 P12) (A. 14) 

where tr^ is the operator tr^ in the space 1. Equation flA.141) just shows that the partial 
trace is not cyclic. As a consequence, one cannot prove that transfer matrices with different 
spectral parameters commute, and the integrability of the model is not guaranteed. 



B Twisted version of XX and Hubbard models 

The Hubbard models we have constructed depend on a single free parameter U . We present 
here a construction that allows us to introduce more parameters. In particular, we will obtain 
an Hermitian Hamiltonian that depends on phases that can be identified with a Aharonov- 
Bohm phase. Again, we present in detail the construction for XX models, and just sketch 
the generalization to Hubbard models. 

We start with an universal XX model defined on the vector space V, with a projector 7i 
such that 7r(V) = W. 

Definition B.l To any projection tc, a refinement is a decomposition 

7T = ©J =1 7r (i) with 7r (i) 7r (/) = 5 j:f vr (i) V j, f = 1, . . . , d 
For each refinement of n and 7f, we introduce 

d d 

JT 12 = I ® I + J2 J2(q a a - 1) vr (a) <g> 

a=l a=l 

where q a& are some non-zero complex numbers. Note that q a5i 7^ ensures that JF 12 is 
invertible. Its inverse reads 



3 In fact it is true for any operator O obeying O(Tn) C J-n- 
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The twisted version of the XX model is defined by the R-matrix 

R12 (A) = T\ 2 #12 (A) T21 
where i?i2(A) is the matrix (12. 4p . We remind that it depends on the choice of the projector 



TV. 



Property B.2 The R-matrix Ri 2 (X) can be rewritten as 

Ru(X) = S12 sin A + (E12 + (I <8> I - S12) cos a) P 12 (B.l) 

Si 2 = ^i 2 Si2^ 21 1 = EEf?a^ (a) ®^ ) + — Tr^OTrW) (B.2) 

1 - 1 V 9ia / 

a=l a=l x ' 

It obeys all the properties stated in theorem \2.1\ but the symmetry one. 

Proof: The expressions (IB. 1[) and (IB. 2ft follow from direct calculations. 
Remarking that for any diagonal matrix P12, we have the property 

^13^23 #12 (A) =R l2 (A) ©13^23 

the YBE for i? 12 (A) is deduced from the YBE for i? 12 (A). For instance, starting from the 
l.h.s. and using the notation \j = \ — \f 

-Rl2(Ai2) -Rl3(Al3) #23^23) = •^12 #12^12) F\ 21 F 13 #13^13) T 31 T 23 #23^23) F 32 
= T\2 i?12(Al2) TnT^^T T 2 \ Rvi^iz) F n -^23 #23(^23) F 32 1 
= T\i Rni^n) -Rl3(Al3) T 31 J? 7 2 \ T 23 #23^23) 32 1 

= T\2 !F\z T-2Z Rvii^Vl) -Rl3(Ai3) #23^23) F\ 31 F 21 

The r.h.s. is treated in the same way. 
Regularity is obtained as follows 

#12(0) = T\2 #12(0) T 2 \ = T 12 P\2 F 2 i = T\ 2 T l2 P\ 2 = P\ 2 
Similar calculations lead to the other properties. ■ 

Property IB .21 ensures that models based on twisted /^-matrices are also integrable and 
possess a local Hamiltonian. Indeed, the corresponding Hamiltonian reads 

L 

H(q a a) =^2Hj,j+i(qaa) with H jd+1 = P jtj+1 E jd+1 (q a - a ) (B.3) 
i=i 

Its hermiticity depends on the parameters q aa - From the calculation 

L 

H(q aa ) = ^2Pjj +1 %j +1 (qaa) (B.4) 
L L 

H(qaa) ] = ^E iJ+ i(g^)P^ + i = ^P iii+1 E i+1 j(0 (B.5) 
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and the identity 

one deduces that H is hermitian when the parameters are phases: 

Qaa = e 4 ^ , 9 aa eR. 

Chosing the parameters q aa to be phases, we have in this way a general Hermitian Hamil- 
tonian with Aharonov-Bohm phases on each site of the model. 

Twisted Hubbard models: The same construction can be done for Hubbard models, 
with now 

Rl l 2 (X)=^U(q\q l )Rl l 2 (X) (^teW))" 1 with ^(q\q l )=^ 12 (q r )^ 2 (q l ). 
Since T and C are diagonal, the new i?-matrix reads 

flg(Ai, A 2 ) = R\ 2 (X 12 ) R{ 2 (\ 12 ) + ^4 tenh^y i?l 2 (A' 12 ) Cj %(A' 12 ) C{ (B.6) 
where R £ (X), e =f , I, have the form (IB. II) . It leads to an Hamiltonian 

L 

that is hermitian as soon as the parameters q e an are phases. 
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